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Abstract 

This contribution is dedicated to dilucidating the role of the Gibbs entropy 
in the discussion of the emergence of irreversibility in the macroscopic world 
from the microscopic level. By using an extension of the Onsager theory to 
the phase space we obtain a generalization of the Liouville equation describing 
the evolution of the distribution vector in the form of a master equation. This 
formalism leads in a natural way to the breaking of the BBGKY hierarchy. 
As a particular case we derive the Boltzmann equation. 



1 



Typeset using REVTgX 



PACS numbers: 05.20.-y, 05.20.Dd, 05.70.Ln 

Introduction. — According to the mechanicistic interpretation of the physical world, the 
basic laws of nature are deterministic and time reversible. However, at the macroscopic 
level we observe irreversible processes related to energy degradation, which generate entropy. 
How do we reconcile the 'spontaneous production of entropy' with the time reversibility of 
the microscopic equations of motion?. At the end of nineteenth century Boltzmann tried 
to answer this question from a probabilistic point of view. According to him, entropy is 
a measure of the lack of knowledge of the precise state of matter, and can be defined as a 
function of the probability of a given state of matter, i. e. it is a function of the microstate. All 
systems in their irreversible evolution tend to a state of maximum probability or maximum 
entropy -the state of equilibrium-. 

In contrast to the Boltzmann entropy, the Gibbs entropy is not a function of the indi- 
vidual microstate but a function of the probability distribution in a statistical ensemble of 
systems, both coinciding at equilibrium. As a consequence of the incompressible character 
of the flow of points representing the natural evolution of the statistical ensemble in phase 
space , the Gibbs entropy is a constant of motion. Thus, it has been argued that the relevant 
entropy for understanding thermodynamic irreversibility is the Boltzmann entropy and not 
the Gibbs entropy [1], [2], [3]. 

But in our opinion, this is not the end of the story. It is our contention here to show 
that depending on the representation we use to describe the state of the system, the Gibbs 
entropy is a good definition of the nonequilibrium entropy which increases in the approach 
to equilibrium and is compatible with the Boltzmann's account for irreversibility . To 
accomplish that goal, we will use the methods of mesoscopic nonequilibrium thermodynamics 
MNET [4], [5] in addition to the generalized Liouville description of dynamical systems in 
terms of the distribution vector. 

The paper is organized as follows. In Section II we introduce the representation of the 
state of the isolated system in terms of the hierarchy of reduced distribution functions. Sec- 
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tion III, is devoted to developing the thermodynamic analysis, and to deriving the entropy 
production. Here, we draw kinetic equations as one of the consequences, in particular the 
Boltzmann equation. In section IV, we stress our main conclusions. 

Distribution vector dynamics. — Let us think of a dynamical system having N interacting 

degrees of freedom (q 1 , ,Q2,Pi, ,Pn), where and pi are the generalized coordinate 

and conjugated momentum corresponding to the i-th degree of freedom. Let H({q N ,p N }) 
be the Hamiltonian of the system, given by 

N ( 2 1 1 N 

^ = E{^ + ^)} + 2 E/(I^-^D > (!) 

where m is the mass, 4> e (qj) = 4> e j is an external potential, and 4>(\qj — qk\) = 4>jk is the 
interaction potential. The state of the system is completely specified at a given time by 
the N-particle distribution function F ({q N ,p N } ; tj, which evolves in time according to the 
Liouville equation. Nonetheless, an alternative description of the state of the system can be 
given in terms of the distribution vector [6] 

i = {foJi(x 1 )J 2 (x 1 ,x 2 ), Jn(xi,x 2 , ,x N )} , (2) 

which represents the set of all the reduced s-particle distribution functions fs(xi,x 2 , , x s ) 

(s = 0, , N), defined through 

N\ r 

f s = jjy _ s y J F(xi, ,x s , x s+ i, , x N ) dx s+ i dx N , (3) 

where Xj = (qj,pj). This function is normalized according to 

r W 

J fsdx, dx s = (jV _ g) , • (4) 

The evolution equations of the s-particle distribution functions can be obtained by integrat- 
ing the Liouville equation constituting a set of coupled equations: the BBGKY hierarchy. 
An arbitrary component equation of this hierarchy is 
d ( s s \ 

7rJs{xi, ,X S )=\Y,L°+ E L 'j,n ] fs(xi X a ) 



dt 



3 ' ^ Jy 

\j =1 j<n=l 



+ E / L 'j,s+Js+i(xi, -;X S+1 ) dx s+1 , (5) 



which shows that the time derivative of f s depends both on f s and on the higher-order 
functions f s +i- Here 



I'l- H°,... 



v 

where [..., ...] p is the Poisson bracket, H° = ^- + <f>p and 



(6) 



J I TTl 



P 



(7) 



with H' j n = \4> jin . 

From Eq. (5) one can infer the evolution equation for the distribution vector 

j t m = cm , 



(8) 



which constitutes the generalized Liouville equation, and £ is the generalized Liouvillian. 
This equation should be equivalent to Eq. (5), thus, the matrix elements {s \C\ s') must have 
the following values: 



(s\C\s') = V<|E^°+ E 

j=l j<n=l 



+S s ', s +i J f^2L' j s+1 ^dx s+1 , 



(9) 



and 



{0\C\s'} = , 



(10) 



where | s) represents the s-particle state. At this point, it is convenient to introduce the 
operator V providing the diagonal part of a matrix , and Q, its complement which gives the 
nondiagonal part. Thus, by using these operators we can rewrite Eq. (8) as follows 



j t - VC\ f it) = Q£S(t) , 



(11) 



with 



i j=l j<n=l 



(12) 



and 



(s \QC\ s') = 5 s > jS+1 




(13) 



The left-hand side of Eq. (11) is invariant under time reversal whereas this symmetry is 
broken in the right-hand side term, as one can convince oneself changing t by —t and taking 
into account Eqs. (12) and (13). Thus, the diagonal part of the generalized Liouvillian 
accounts for the reversible evolution of the distribution vector, and the nondiagonal part 
introduces irreversibility in the dynamics of f leading to dissipation, this fact will be show 
adequately in the section that follows. Hence, in the sense that the time reversal invariance 
is broken, the generalized Liouville equation given through Eqns. (8) or (11) is irreversible. 
Irreversibility is manifested in the dynamics of the system when we use the adequate de- 
scription, i.e. in terms of the distribution vector, this description brings with itself certain 
degree of coarse graining. 

Mesoscopic nonequilibrium thermodynamics. — The MNET constitutes the mesoscopic 
generalization of the macroscopic formalism of the non-equilibrium thermodynamics [8]. 
Here, the entropy is a functional of the distribution vector in the phase space given by the 
Gibbs entropy postulate [7] . Moreover, we apply the principles of equilibrium thermodynam- 
ics at a local level in phase space, and we follow the general scheme of non-equilibrium ther- 
modynamics in order to compute the entropy production from the Gibbs entropy postulate. 
One then postulates the linear phenomenological equations relating fluxes and thermody- 
namic forces occurring therein. Immediately afterwards, by substituting this phenomenolog- 
ical equations in the generalized Liouville equation written in the form given in Eq. (11), one 
derives master equations for the distribution vector. Therefore, the entropy of the system, 
in analogy to the Gibbs entropy postulate, is given by 



where ks is the Boltzmann constant, S Q the equilibrium entropy, and f Q is assumed to 



S = -k B Tr{Uog (f^f)} + S 




(14) 
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be the equilibrium distribution vector satisfying £f G = 0. The rate of change of S can be 
obtained by differentiating Eq. (14) 

which for small deviations from the equilibrium reduces to 

£ = - kBTr (16) 

Here, the only contribution to the entropy production comes from the irreversible change of 
the distribution vector since the change reversible contributes to zero entropy production. 
The reversible change of the distribution vector is given by the mainstream current induced 
by the s-Hamiltonian flow in each one of the s-particle subspaces of the total phase space, 
and it is known that the Hamiltonian flow does not create entropy. According to the Onsager 
theory [8], [9] from Eq. (16) we infer the linear relation 

- = M{f- 1 (f c -f)}=MX, (17) 

where the master or mobility matrix M acts on an arbitrary vector Y according to 

(s|M|s')(s' | Y =j M(s\s')Y s/ (x' 1 ,....,x' s ,)dx' 1 dx' s , . (18) 

Substituting the linear relation Eq. (17) into Eq. (15) and using the cyclic invariance of the 
trace, we find for the entropy production 

Tr(XMX)>0, (19) 

according to the second law of Thermodynamics. This relation imposes the Hermitian 
character of the master matrix 

M(s\s') = M(s'\ S y , (20) 

as predicted by the Onsager symmetry relations. Here ' refers to the Hermitian conjugate. 
Furthermore, due to the fact that f should be normalized Tr(f) is a constant quantity which 
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is a function of N, from Eq. (17) we infer that Tr(MX) = 0, which, when we take into 
account the Hermitian character of the master matrix, leads to the following constraints 

J M(s | s')dx 1 dx s = J M(s | s'fdx[ dx' s , = . (21) 

Hence, by using the generalized Liouville equation (8) and Eq. (17) it comes to be that 
Eq. (11) can be written in a more convenient form 

(!- p£ ) f w- s fH Mx - < 22 > 

with 

(s\QM\s')(s' | X 

= 8s> t s+\ j M(s + l\s')X s ,dx' 1 dx' s ,dx s+1 , (23) 

written in analogy with Eq. (13). It should be mentioned that the presence of the master 
matrix in Eq. (22) notably simplifies the BBGKY hierarchy. In fact, the master matrix 
introduces a relaxation time scale. Hence, our theory is equivalent to a relaxation time 
approach to the study of the BBGKY hierarchy valid when there is a broad separation 
between the hydrodynamic and microscopic scales. 

Here, Eq. (22) can be brought into a more common form, by introducing a new function 
W defined through [10] 

fo,s>W{s'\s) =-M(s\s') 

+5(x[-x 1 ) S(x' s , - x s )tp s (x 1 , ....,x s ) , (24) 

where the auxiliary function ij) s is not arbitrarily selectable because of the constraints given 
by Eq. (21). Rather applies 

*Ps (xi, , x s ) = f ,s(xi, , x a ) J W(s\s')dx[ dx' s , . (25) 

Thus, Eq. (22) may be written in terms of the transition matrix as follows 

- VC\ f (t) = QWf (t) , (26) 



where 



{s\QW\s')(s'\f=S al , a+1 J{fAx\ 



,x' s ,)W(s'\s + l) 



x s+1 )W(s + 1 | s')}dx[ dx' s ,dx. 



s+l ■ 



(27) 



It should be stressed that the introduction of the transition matrix leads to a deeper physical 
insight of our result. 

A particularly interesting case corresponds to s — 1 where Eqs. (26) and (27) reduce to 



This is not yet a kinetic equation; however, its uncorrelated part is and can be written as 



which is the famous Boltzmann equation. The correlated part of Eq. (28) does not matter 
because the correlated part of each side of this equation is identically zero. Analogous 
procedure can be followed to obtain the kinetic equations for the correlations coming from 
the components of f of an order higher than one. 

The quantum generalization of our classical treatment can be done in terms of the Wigner 
distribution vector [6], defined as the set of all the s-particle Wigner functions. This quantum 
approach follows the same lines as the previous classical approach, but with the unavoidable 
differences in the algebraic manipulations. 

Discussion. — Here, we have presented a representation of the statistical description of a 
many-body system in terms of the distribution vector, which reveals the irreversible com- 
ponent of the motion of this system. This component is hidden when we work with the 
standard Liouville equation but at the level of the generalized Liouville equation which en- 
tails a certain degree of coarse graining, this irreversible component becomes obvious. The 
Liouville equation is a closed equation for the phase space distribution function while the 




(28) 




(29) 
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generalized Liouville equation encloses a set of coupled equations, each one in itself repre- 
sents a contraction of the statistical description. But, as we have mentioned in section II, 
the full statistical description in terms of the phase space distribution function is equivalent 
to the description in terms of the set of all the reduced s-particle distribution functions. In 
this scenario, we show that the Gibbs entropy given by the Gibbs entropy postulate is a 
non-conserved quantity throughout the motion in phase space. This entropy should increase 
until its maximum value at equilibrium according to the second law. The consequences we 
have drawn from the Gibbs entropy are compatible with the explanation of irreversibility 
which follows from Boltzmann's entropy postulate. 

By applying Onsager's theory in phase space, i.e. in the framework of MNET, we obtain 
a linear relation between the irreversible rate of change of the distribution vector and the 
conjugated thermodynamic force, which introduces a master matrix which should satisfy 
the properties required by the Onsager theory. This result simplifies the BBGKY hierarchy 
and enables us to break the hierarchy obtaining the kinetic equation for the one-particle 
distribution function which coincides with the Boltzmann kinetic equation. Hence, this last 
result constitutes a test of our theory. Hence, we have shown a method based on the MNET 
for deriving kinetic equations. 

We also stress that in our mesoscopic approach, the master matrix introduces a relaxation 
time scale which shows that our theory is valid when there is a broad separation between 
the hydrodynamic and microscopic scales. This separation of time scales allows the rate of 
change of the distribution vector due to the main stream flow velocity in phase space balance 
the irreversible change in order to achieve a stationary distribution at long time. This is 
precisely what happens in the context of the phenomena described by the kinetic equations. 

I want acknowledge Prof. J.M. Rubi for profitable discussions. This work was supported 
by DGICYT of the Spanish Governement under Grant No. PB2002-01267. 
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